function [x,k]=Jacobimethod(~,~,~,~,~)%A是线性方程组的左端矩阵
%b是右端向量
%x0是迭代初始值
%N表选代次数上限，若迭代次数大于N，则选代失败% emg 表示控制精度
%用Jacobi迭代法求线性方程组Axb的解
%k表迭代次数
%x表示用选代法求得的线性方程组的近似解n=length(A);
xl=zeros(n,1);x2=zeros(n,l);
x1=x0;k=0;
r=max(abs(b-A*xl));
while r>emg
for i=l:n
sum=0;
for j=1:n
if i~=j
sum=sum+A(i,j)*x1(j);
end
end